J>1 



(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(19) World Intellectual Property Organization 
International Bureau 

(43) International Publication Date 
21 February 2002 (21.02.2002) 




PCT 



(10) International Publication Number 

WO 02/14934 Al 



(51) International Patent Classification 7 : G02F 1/01 

(21) International Application Number: PCT/US0 1/1 9788 

(22) International Filing Date: 21 June 2001 (21.06.2001) 



(25) Filing Language: 

(26) Publication Language: 



English 



English 



(30) Priority Data: 
00401993.1 



1 1 August 2000 (1 1 .08.2000) EP 



(71) Applicant: CORNING INCORPORATED [US/US]; 1 
Riverfront Plaza, Corning, NY 14831 (US). 

(72) Inventors: FONDEUR, Barthelemy; 6, rue Sergent-Per- 
rier, F-77300 Fontainebleau (FR). RANALLI, Eliseo, R.; 
1406 Reggio Aisle, Irvine, CA 92696 (US). WIG LEY, 
Peter, G.; 2802 Pwoderhouse Road, Corning, NY 14830 
(US). 



(74) Agent: REDMAN, Mary, Y.; Coming Incorporated, SP 
TI 3 1, Corning, NY 14831 (US). 

(81) Designated States (national): AE, AG, AL, AM, AT, AU, 
AZ, BA, BB, BG, BR, BY, BZ, CA, CH, CN, CR, CU, CZ, 
DE, DK, DM, DZ, EE, ES, FI, GB, GD, GE, GH, GM, HR, 
HU, ID, IL, IN, IS, JP, KE, KG, KP, KR, KZ, LC, LK, LR, 
LS, LT, LU, LV, MA, MD, MG, MK, MN, MW, MX, MZ, 
NO, NZ, PL, PT, RO, RU, SD, SE, SG, SI, SK, SL, TJ, TM, 
TR, TT, TZ, UA, UG, UZ, VN, YU, ZA, ZW. 

(84) Designated States (regional): European patent (AT, BE, 
CH, CY, DE, DK, ES, FI, FR, GB, GR, IE, IT, LU, MC, 
NL, PT, SE, TR). 

Published: 

— with international search report 

For two-letter codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begin- 
ning of each regular issue of the PCT Gazette. 



(54) Title: TAPPED DELAY LINE BASED GAIN FLATTENING FILTER 

44 



m 

ON 

O 




(57) Abstract: A method of modifying the spectral distribution of an optical signal includes the steps of splitting the signal among 
a plurality of optical paths (14), delaying the propagation of the signal through one or more of the paths by a time duration which 
results in the desired spectral distribution when the signals at the output ends of each optical path are recombined, and recombining 
those outputs. An apparatus for changing the spectral profile of an optical signal includes an input waveguide and a plurality of 
tapped delay lines (24, 26, 28, 30). Each tapped delay lines has a delay element (38) configurable to impart a predetermined delay to 
the signal propagating down the delay line. Couplers (18, 19, 20) split an optical signal, propagating on the input waveguide, among 
the tapped delay lines (24, 26, 28, 30). The delays imparted to the signals on each line are chosen so as to effect the desired change 
to the spectral profile. 
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TAPPED DELAY LINE BASED GAIN FLATTENING FILTER 
Background of the Invention 

1. Field of the Invention 

The present invention relates generally to optical filters, and more 
specifically to configurable filters for spectral shaping. 

2. Technical Background 

Actively configurable optical filters have many potential applications in optical 
communications, such as temporal pulse compression, dispersion slope compensation 
and spectral shaping. Spectral shaping is of particular interest in connection with 
optical amplifiers. The gain profile of a typical optical amplifier may be affected by 
numerous factors such as, for example, fiber length, fiber composition, splice losses, as 
well as the population inversion level of the dopant atoms (e.g. erbium) providing the 
gain. Having a relatively large optical signal at the input of the amplifier tends to 
deplete the excited electron population, resulting in a different spectral gain profile than 
that for a small input signal. Average gain also varies with input signal power. 

One approach to dealing with this effect is to operate the amplifier at a constant 
inversion population, with a fixed spectral filter that provides gain flattening over the 
operating wavelength range. In this approach, the inversion level is held fixed by 
selectively attenuating the input signal until the desired inversion level is achieved. 
However, there may be circumstances where operating the amplifier at a constant gain 
level, rather than a constant inversion level, would be advantageous. By utilizing an 
actively configurable spectral filter, an amplifier could be used in a constant-gain mode, 
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rather than constant inversion level, by flattening the gain profile of the amplifier for 
arbitrary inversion levels. One approach for realizing such a configurable spectral filter 
is to spatially disperse the various wavelength components of the input signal via a 
diffraction grating, then spatially modulate the dispersed components via a spatial light 

5 modulator or acoustically-formed phase grating. The polarization dependence and 
diffractive losses of the grating could be a drawback in some applications, however. 
Prisms could be used instead of a diffraction gratihg'to eliminate polarization 
dependence and to decrease angular dispersion. Prisms require an increase in package 
size, however, which makes their use undesirable in many applications. 

10 What is needed, then, is an actively configurable spectral filter that is 

polarization-independent, which can fit in a relatively small package, and which does 
not suffer from diffractive losses. 

Summary of the Invention 

15 One aspect of the invention is an apparatus for changing the spectral profile of 

an optical signal that includes a plurality of optical paths, each having associated with it 
a delay element which imparts a predetermined delay to a signal propagating through 
the optical path. A plurality of couplers is configured to split the optical signal among 
the plurality of optical paths according to a predetermined splitting ratio. The 

20 predetermined delays and predetermined splitting ratios are chosen so as to effect the 
desired change to the spectral profile. 

In another aspect, the delay elements may include portions of the waveguide 
which have an altered index of refraction. The index of refraction may, in yet another 
aspect, be reversibly alterable. 

25 In another aspect, the invention includes a method for modifying the spectral 

distribution of an optical signal which includes the steps of splitting the signal among a 
plurality of optical paths, each having an output end, so that a portion of the signal 
propagates on each of the optical paths; changing the relative phases of the portions of 
the signal propagating through the optical paths; modifying the amplitudes of the 

30 portions of the signal propagating through the optical paths; and combining the portions 
of the signal at the output ends of the optical paths. The relative phase changes and the 
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amplitude modifications have magnitudes that result in a desired spectral distribution 
when the portions of the signal at the output ends of the optical paths are combined. 

In another aspect, the invention includes a method for modifying the spectral 
distribution of an optical signal which includes the steps of splitting the signal among a 

5 plurality of optical paths according to predetermined splitting ratios, and changing the 

relative phases of the portions of the signal propagating through different ones of the 
optical paths; the relative phase changes and the splitting ratios having magnitudes 
which result in a desired spectral distribution when the portions of the signal at the 
output ends of the optical paths are combined; and combining the portions of the signal 

10 at the output ends of the optical paths. 

In another aspect, the invention includes an optical communication system 
comprising an optical amplifier and a spectral filter that includes an input waveguide 
and a plurality of tapped delay lines. Each of the tapped delay lines has a delay element 
configurable to impart a predetermined delay to the signal propagating down the delay 

15 line. Couplers split an optical signal, propagating on the input waveguide, among the 

tapped delay lines. The delays imparted to the signals on each line are chosen so as to 
effect the desired change to the spectral profile. 

The claimed embodiments of the invention provide spectral shaping which is 
polarization-independent, and which, in certain embodiments, may be actively 

20 configurable or permanently configured. The spectral filtering of the invention may, if 

desired, be provided in a planar device, which is compact, reliable and economically 
fabricated. It is particularly well suited to use as a gain-flattening filter for an optical 
amplifier, which allows the amplifier to be used in a constant gain mode. 

Additional features and advantages of the invention will be set forth in the 

25 detailed description which follows, and in part will be readily apparent to those skilled 

in the art from that description or recognized by practicing the invention as described 
herein, including the detailed description which follows, the claims, as well as the 
appended drawings. 

It is to be understood that both the foregoing general description and the 

30 following detailed description are merely exemplary of the invention, and are intended 

to provide an overview or framework for understanding the nature and character of the 
invention as it is claimed. The accompanying drawings are included to provide a 
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further understanding of the invention, and are incorporated in and constitute a part of 
this specification. The drawings illustrate various embodiments of the invention, and 
together with the description serve to explain the principles and operation of the 
invention. Additional features and advantages of the invention will be set forth in the 

5 detailed description which follows, and in part will be readily apparent to those skilled 
in the art from that description or recognized by practicing the invention as described 
herein, including the detailed description which follows, the claims, as well as the 
appended drawings. 

It is to be understood that both the foregoing general description and the 

10 following detailed description are merely exemplary of the invention, and are intended 
to provide an overview or framework for understanding the nature and character of the 
invention as it is claimed. The accompanying drawings are included to provide a 
further understanding of the invention, and are incorporated in and constitute a part of 
this specification. The drawings illustrate various embodiments of the invention, and 

15 together with the description serve to explain the principles and operation of the 
invention. 



Brief Description of the Drawings 

Figure 1 is a schematic diagram of an embodiment of the invention; 
20 Figure 2 is a schematic diagram of an alternate embodiment of the invention; 

Figure 3 is a schematic diagram of another alternate embodiment of the 
invention; 

Figure 4 is a schematic diagram of another .alternate embodiment of the 
invention; 

25 Figure 5 is a graph showing a gain spectrum for an EDFA before gain 

flattening; and 

Figure 6 is a graph showing a gain spectrum for an EDFA after gain flatteneing. 

Detailed Description of the Preferred embodiments 
30 Reference will now be made in detail to the present preferred embodiments of 

the invention, examples of which are illustrated in the accompanying drawings. 
Wherever possible, the same reference numbers will be used throughout the drawings 
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to refer to the same or like parts. An exemplary embodiment of an apparatus for 
changing the spectral profile of an optical system of the present invention is shown in 
Figure 1, and is designated generally throughout by reference numeral 10. 

In accordance with the invention, an embodiment of the invention provides 
5 spectral shaping by employing an optical tapped delay line realization of the desired 

narrow-band impulse response. In this embodiment, a splitter tree partitions the 
incoming signal among several taps in the correct proportions for a given impulse 
response. Each tap imparts a nominal delay of T seconds, relative to neighboring taps. 
The design of the splitter tree may be chosen so that the inverse of the delay (1/T) is on 

10 the order of the free spectral range of interest, or more preferably, about twice the 

desired free spectral range. Delay perturbations in each tap set the approximate phase 
for that tap. With the splitter ratio determining amplitude, and the delay perturbations 
determining phase for the signal propagating down each tap, the signals output from 
each tap are superimposed to yield an output signal having the desired spectral shape. 

15 As illustrated in Figure 1 , an embodiment of the present invention for changing 

the spectral profile of an optical signal includes an array of optical paths 14 which, in 
the illustrated embodiment, are configured in a splitter tree arrangement. The optical 
path array 14 can be an array of waveguides in a planar substrate, such as, for example, 
single mode silica waveguides on a silicon substrate, or waveguides on a LiNi03 

20 substrate. In the illustrated splitter tree arrangement, a first optical path 16 directs the 

signal to a first coupler 1 8, which may be, for example, a Mach-Zehnder 
interferometer or other suitable variable coupler so that the coupling ratio may be 
selectively configured as desired. The first coupler .1 8 splits the signal into two 
portions and couples the two portions into optical paths 21 , 22, each of which is 

25 optically connected to one of a second and a third coupler 19, 20 respectively, which 

may be any suitable variable 2X 2 coupler such as a Mach-Zehnder interferometer. 
Each of the second and third couplers 19, 20 split the power input to it on to a pair of 
waveguides 24, 26 and 28, 30, respectively. In the illustrated embodiment, these 
waveguides act as taps, with the relative tap weights being effectively set by splitting 

30 the input power in appropriate ratios over the two stages of variable couplers made up 
by the first coupler 1 8 (in the first stage) and the second and third couplers (in the 
second stage). Each of the tap waveguides 24, 26, 28 and 30 has a different effective 
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optical path length chosen so that each tap imports a nominal delay T to a signal 
propagating through it, relative to a neighboring tap. As illustrated, this may be 
accomplished by the tap waveguides having successively longer (or shorter) actual 
lengths. For purposes of illustration, a two-stage splitter tree is illustrated in Figure 1. 
5 Embodiments with different numbers of stages and of tap waveguides are within the 
scope of the invention. 

A delay element 38 is associated with each of the tap waveguides 24, 26, 28 and 
30 for imparting delay perturbations to the Signal traveling through that tap. The delay 
elements 38 may be, for example, portions of the associated waveguide which have an 
1 0 altered index of refraction. Altering the index of refraction in a waveguide changes the 
effective optical path length through that waveguide, thereby effecting a delay in 
propagation time, and in phase, through that waveguide. The alteration in refractive 
index can be a permanent alteration, or alternatively, a temporary or reversible 
alteration. In this alternative, the refractive index preferably could be selectively 
1 5 modified as desired to achieve the desired delay in propagation. 

Some examples of how this modification of refractive index can be achieved 
include, for waveguides in a silica substrate, exposing a portion of the waveguide to 
UV radiation. This will result in a permanent modification of the refractive index. 
Temporary or reversible modifications to the refractive index can be provided for 
20 waveguides in a silica substrate by providing selectively actuated heating elements 

adjacent the waveguide portion for which the refractive index changes are desired. For 
an embodiment of the invention using waveguides in an electro-optic substrate such as 
LMO3, electrodes adjacent a waveguide may be used to provide a voltage across the 
waveguide, which will change the refractive index in a non-permanent fashion 
25 according to the well-known electro-optic effect. 

In the illustrated embodiment, the tap waveguides are reflected back on 
themselves by, for example, a reflector 34 at the output end of the optical paths which 
make up tap waveguides 24, 26, 28 and 30, reflect the portions of the signal 
propagating down these optical paths back through the splitter tree 14 in the opposite 
30 direction. This ensures that the couplers, which separate the input signals into the 

successive stages of the splitter tree and the coupling of the delayed copies of the input 
signals, are accomplished in the same ratios. The reflector 34 can be a polished gold- 
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plated surface or any other suitable reflective optical element or elements. The ports 
40, 42 and 44 to the first, second and third couplers 18,19 and 20, which do receive 
input power may be used to probe the couplers, for system monitoring. 

Figure 2 shows an embodiment of the invention similar to that of Figure 1, in 

5 which the first, second and third couplers 1 8, 19, 20 are illustrated as Mach-Zehnder 

interferometers 41, 43, 45, respectively. In an embodiment where the splitter tree 14 is 
fabricated with silica waveguides on a silicon substrate, the splitting ratios through 
these Mach-Zehnder elements 41, 43, 45 may be modified as desired by actuation of 
heating elements 46, 47, 48 adjacent the waveguide making up one arm of the Mach- 

10 Zehnder interferometer. In the embodiment of Figure 2, a circulator 50 between the 
input waveguide 52 and the splitter tree 14 ensures that the reflected signal propagates 
to the output optical path 54. Those skilled in the art will appreciate that alternative 
means may be used to direct optical signals from a communication system into and out 
of the device 10, employing other optical components such as isolators, for example. 

1 5 The embodiments of Figures 1 and 2 may be modified by replacing the reflector 

element 34 with a mirror image of the array. An example of this is illustrated in Figure 
3. 

In another aspect of the invention, illustrated in Figure 4, an optical 
communication system includes a transmitter 53, a receiver 54, an optical amplifier 56 

20 which may be an erbium-doped fiber amplifier (EDFA), and a gain-flattening filter 58. 

The gain-flattening filter can be a spectral filter 1 0 as illustrated in Figures 1 , 2 or 3. A 
waveguide 60 propagates a signal from the amplifier 56 to the gain flattening filter 58. 
In this aspect, the splitting ratios and delays for the embodiment illustrated in Figure 1, 
2 or 3 would preferably be chosen to provide a gain flattening effect to the output of the 

25 optical amplifier 58. 

For example, suppose it is desired to realize a gain-flattening filter for an 
erbium-doped fiber amplifier (EDFA), operating over the C-band (roughly 1530-1560 
nm), with 8 taps. In this particular example, the parameters will be chosen so that the 
device optimally flattens the gain of the amplifier at a population inversion level of 

30 65%, the expected operating condition for a particular amplifier design. Figure 5 

shows the spectral gain profile for the amplifier at that inversion level. It is recognized 
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that the spectral measurement is noisy, and the optimization should be robust enough to 
accommodate this limitation. 

Calculation of the appropriate complex tap weights may be done using 
MathCAD software, and an annotated MathCAD script is given in Appendix - Part 1 to 
5 illustrate an example of appropriate steps taken in designing a device to realize such a 
function. Spectral scans at several inversion levels are recorded in the file, C-Band.txt, 
from which the design wavelengths (indicated bythe'vector k - Design) and gains 
(indicated by the matrices: y in dB units, Ain absolute units) are extracted. 

The Nyquist (sampling) theorem ordinarily dictates that the sampling interval T 

10 be chosen to correspond to the inverse of twice the frequency range of interest. In the 
present design example, this criterion would lead to the specification T = 1 10 fsec. 
However, this should not be interpreted as a rigid constraint, and in fact a minimum 
least square error solution is actually obtained from 185 fsec. In fact, if the Nyquist 
theorem is followed blindly, there is also the very real possibility that a quite lossy 

15 system will result, since the range over which the desired response is specified is only 
half the free spectral range (in fact, it is this condition which insures that there is 
absolutely no aliasing), and over the other half may be where the transmission 
maximum occurs. This possibility follows from the constraint that a given tap weight 
cannot exceed unity. Thus, some searching for an optimal T must occur in general. For 

20 a given T however, optimization proceeds as follows. 

The first step is to quickly arrive at a linear least mean-square (LMS) solution. 
This leads to a solution for the complex tap coefficients which is obtained by 
representing the system transfer function as a weighted (by the tap coefficients) 
superposition of sinusoids in the frequency domain. Under such an optimization 

25 paradigm, the transfer function to be approximated must be specified both in amplitude 
and phase. Since the design task simply specifies the amplitude of the desired spectral 
response, the phase can be any arbitrary function. The phase functions <D in Appendix - 
Part 1 were chosen based on a minimum-phase criterion. In the same way that the real 
and imaginary parts of the refractive index are related in dielectric materials, the 

30 principle of causality leads to this relationship between the amplitude and phase of a 

linear, causal system. Clearly, this is not the only phase relationship that can exist for a 
given amplitude of response, and in fact a non-linear optimization leads to a different 
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function, as will be explained shortly. However, it is necessary to initiate the 
optimization at a reasonable starting point, and in this example, this function provides 
that initial guess. 

Thus, with amplitude and phase determined, the target transfer function H is 
completely specified. The graphs of X_ design k versus |H k, 3I and <J> k,3 illustrate the 
amplitude and phase for this function, chosen for N=3, where the population inversion 
is 65 %. In the present example, the number of taps, indicated by the constant, Taps, is 
8. The complex cumulants are then computed for the LMS estimation of the complex 
tap weights, represented by the vector h, which is then normalized and presented 
graphically in the MathCAD script. From this, two frequency responses are derived: 
H2 and H3. H3 represents the response of an ideal tapped delay line system, for which 
h is not frequency dependent, while H2 represents the response for a system in which 
the phase is allowed to vary over the frequency range in a way consistent with 
approximating a desired phase by a constant delay over a narrow range of frequencies. 
The graphs of AJData m verses |H 2 ml and |H 3 ml show the difference between the two 
functions. There is little difference between the functions, validating the assumption 
that for this narrow-band system, a constant delay is a good approximation to a constant 
phase. Over the frequency range of interest in the example, the RMS ripple in the 
compensated output is only 0.2 dB. 

As mentioned above, the linear prediction may be used as a starting point for 
optimization for a given design. Once h is found in this way, a non-linear optimal 
solution may be found by entering the elements of h into the C program shown in 
Appendix - Part 2, which randomly perturbs the real and imaginary components of h in 
such a way as to minimize the RMS error in the gain-compensated ripple, after an 
appropriate number of trials and with simulated annealing. 

Following this non-linear optimization, the optimal impulse response 
coefficients (i.e. the elements of h), given in Appendix-Part 3, show that the RMS 
ripple has been reduced to a mere 0.047 dB. As can be seen in Figure 6, this ripple is 
most likely due to random error in the data. 

It will be apparent to those skilled in the art that various modifications and 
variations can be made to the present invention without departing from the spirit and 
scope of the invention. Thus, it is intended that the present invention covers the 
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modifications and variations of this invention provided they come within the scope of 
the appended claims and their equivalents. 
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Appendix -Parti 

Data which, in this example, consists of gains for 71 wavelengths at 1 1 inversion levels, is read in: 
Spectra="C-Band.txr 

5 V := READPRN( Spectra) 

M:=71 N:=11 

m := 0~ M - 1 
n :* 0-N - 1 
10 Data wavelengths 



15 



25 



30 



X_Data^ '•- V ^ n 
— m m t i ,u 



Inversion levels 



n2 n :=V 0.n+.15 



Design wavelengths. In this example, only the useful portion of the gain curves are of interest. 

k := 0.34 
20 XJDesign k := X_Data R 

Radian frequencies: 

2-7I-3.10 8 



X Data .10* 9 
— m 



<3ains at kth wavelength, nth inversion level, in dB 

y k.n = V k+I.n + 15 

Corresponding sampling rate, in seconds. 

T:= 1.85-10" 13 
Absolute spectral response. 



Tk.n 
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Mean radian frequency for n'th inversion level. 



XX. 



Cl_ : 



Create a phase profile, based on a minimum-phase criterion. Assume spectrum to repeat after a free 
spectral range 



s =0.. 34 p := -5.. 5 

p L s 

*k.n :=mod K.n»' t - 2 ) 



35 + k- s)=0,0, 



K A s.n)- 
p-35 + k - s 



Desired transfer function for gain compensation at nth inversion level (choice of numerator sets desired 
gain). 

<yTk.n 



0.06 



| H k,3 | 0-04 



0.02 



i 1 r 




j L 



1520 1530 1540 1SSO 15€ 
Design ^ 



♦k,3 O - 




-O.S 



1520 1S30 1540 1550 1560 1 
A,_Design k 



Number of waveguide taps (in this example, 8 taps). 
Taps := 8 
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Choice for an inversion level at which to operate. 

N := 3 
n2 N « 0.65 

Frequency response vector and matrix for least mean square (LMS) error computation, 
t := 0 .. Taps - 1 u := 0 Taps - 1 

k 

k 

Impulse response vector which gives the LMS error, 
h = NT 1 v 

Normalize impulse response. 

Norm:=3 M 
t 

x Norm 
Amplitude t := | h fc | 

2d Amplitude D 1 =1 
t 
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Amplitude and phase of the impulse response which gives the minimum RMS error for the 
chosen number of taps, at the chosen inversion level: 



10 



15 



[h t | 0.5 



h = 



O 5 

0.664294 + 4.84072*10 3 i 
-0.041764+ 0.04456i 
0.107477 0.061672i 
0.0524+ 0.0761131 
-0.018686 + 0.023785i 

-4.532982'10~ 3 - 4.09261 5-1 0~ 3 i 

-4.763722-10"" 3 + 0.012596i 



-7.402329-10 3 " 4.125711*10 ~i 



3. 




20 



Actual frequency response of Tapped Delay Line (TDL) system. 



H2 m := El h tl" exp 
t 



■m 



exp(-jo m -t.T) 



25 



Frequency response of ideal IMS TDL system, 
t 



30 



H2 m | 0.5 



t r-r — T 



i I L 



1520 1540 1560 1580 
J._Data m 



| H3 m | 0-5 



J 1 



1520 1540 1560 1580 160 
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5 Gain compensation: Amplifier + TDL filter: 
Comp k :=10.|og[|A k)N .(]H2 k D 2 |] 



10 



E Comp H 



15 



20 



25 



Normalize impulse response. 

Norm :=X!(l h tD 



h, := — — 
1 Norm 



Ihtl 


arg(hi> 


0.664311 


7.286887-10' 3 


0.061072 


2.323812 


0.123914 


-0.520946 


0.092406 


0.967871 


0.030248 


2.236701 


6.107161 -10" 3 . 


-2.407204 


0.01346 


1.932358 


8.47443 -10" 3 


-2.633122 



H = 24.340649 
max( Comp) - min( Comp) = 0.675054max Ripple in ( 



2(Comp k -ji) 
k = 0.207617 RMS Ripple 



S 1 

k 



30 
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Appendix - Part 2 

/* AGFF_l.c */ 

/* This * program performs non-linear optimization of the complex tap weights' 
for the active gain flattening filter, starting from the LMS predicted 
tap weights from the MathCad simulations. 



ft include <stdio.h> 
ftinclude <stdlib.h> 
ft include <math.h> 



/* Constants */ 



int const Taps = 8; 
double const T • 1.B5E-13; 
int const Trials » 10000; 
double const Step - 0.001; 



/* Sampling interval. 

/* Number of trials in search for optimal tap coefficients. 
/* Step size for coefficient perturbation. 



/* variables */ 

int ra,n,k,LaBt_m; 
float Dummy; 
int Mum; 

int Lambda (1001 ; 

double Omega [100] ; 

double arg[100] ; 

double Gain [100 J ; 

double Target [100] ; 

double Hurdle; 

double h_Ref ; 

FILE * Spectrum ; 

FILE *Outfile; 

double h(16] - { 0.968829, 

_ -0.060909, 

0.156748, 
0.076421, 
-0.027252, 
-0.00661106, 
* -0.0069476, 
-0.010796, 



/* Number of frequencies considered. 

/* Wavelengths in nm. Contained in first column of the file Spec 

/* Corresponding omega* a. 

/* Product of omega and sampling interval. 



/* Successively lowest II*EF values. 

/* Dummy variable for perturbing Re and Xm parts of impulse resp 

/* File contains spectrum to be flattened. 

/* File contains optimized impulse response; 

0.00705988, 

0.064988, 

-0.089945, 

0.111006, 

0.03469, 

-0.005968813, 

0.018371, 

-0. 006017082); 



/* Functions */ 

double S (int Computed&t); 

double XLZFO ; 

mainCargc, argv) 
int argc; 
char **argv; 



{ 



/* Computes the spectrum for the as -then defined 

complex impulse response, at the radian frequency 
indexed by the integer ComputedAt. 

/* Resultant 1LEF. 



Spectrum * f open < "a: \Gain.txt", »r m ) ; 
Outfile - fopent»a.-\h.txf, ? 



/* Read in gain data. First column is wavelength., second is gain, in dB. 

The wavelengths should only cover the range we're interested in flattening. 
After reaching EOF, the variable Hum will represent the number of wavelengths 
represented in this range. 

Nura » 0; 

while (fscanf (Spectrum, "%d %f \n" , RLambda (Num] , &Dummy) la EOF) 

Gain [Num] > pow(io, Dummy /10) ; 

OmegaCNum] « 1. 8 84956E18/ (float) Lambda (Num] ; 
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arg [Nura] » Omega [Hum] *T; 

} 

/* Compute Inverse function for 25 dB overall gain */ 
^ for(m«0;m<Num;nH"h) Target [ml ■ 316/GainCm] ; 

Hurdle - XLEF(); /* initiates search. */ 

Last__m » 0 jm * 0; 

f or (n»0 ;n<Trials ; n++) 

while (m»«Last_m) m - < int) floor (2 - 0 *Taps*rand ( ) / ( float) RAND_MAX) ; 

bjlef m h [ra] # 

him) Step; 

while ( IliBF ( ) <Hurdle) { 

Hurdle m ILEP() ; 
10 h_Ref - htm] % 



15 



20 



}; 

h[m] - hjief ; 
h[m] — Step; 
while (ILEF() <Hurdle) { 



h[m] +- step; 



Hurdle - ILEP() • 
hJRef - h[m) ; 
h[m) — Step; 



html « hjief ; 
Lao t_m - m; 

) 

print f(" %f Nn", Hurdle) ; 

for(m«0;m<Taps;m++) printf(" *f %f \n* ,h[2*m) ,h[2*ra+l] ) ; 
£or(m»0;m<Taps;ra++) fprintf (Out file, » *f tt \n»,ht2*m] ,h[2*m+l] ) ; 

f close (Spectrum) ; 
f close (Outfile) ; 



double S(int At) 

double Dummy_Re, Dummy_Xm; 

DummyJRe m 0 ; "~ . • - • 

Duramy^Xm » 0; 

for (k=0 ;k<Taps ;k++) 

7 * Dummy_Re +» h[2*k) *eos (arg [At] *k) + h[2*x+l] *sin (arg [At] *k) ; 

^ Dummy_Ira h [2*k+l] *cos (arg [At] +k) - h[2*k] *sin(arg(At] *k> ; 

^ return ( Dumniy^Re* Dummy_Jte + Dummy JIm*Dummy_Im) ; 
double ZLEFO 

{ 

double min, max. Dummy ; 
int m; 

min - 10*logl0(S(0)*Gain[0]); 
max » min; 
30 for (m«l;m<Num?m++) 

Dummy *» I0*logl0 (S (m) *Gain [m] ) ; 
if (Dummy?. max) max » Dummy; 
if (Dummy <min) min « Dummy; 

} 

return (max- min.) ; 



WO 02/14934 



18 



PCT/US01/19788 



Appendix - Part 3 * 

Amplitude and phase of the impulse response which gives the minimum RMS error for the 
chosen number of taps, at the chosen inversion level: 



10 



15 



0 5 

0.654353+ 9.581827-10~ 3 i 
-0.03699 + 0.038666I 

0.115147 - 0.071021 

0.054962+ 0.088909i 
-0.020698 + 0.019021 

-2.750918-10" 3 + 1.036431 -lo" 3 i 

~7.045228-10~ 3 + 0.0113451 



-7.764422-1 0* 3 - 1.167252-10 * i 



r3, 



arg o 




20 Actual frequency response of TDL system. 



t 



/ ar 9( h t)\ 
JO- i_S£ ) 



•exp(-j.<D m .t-T) 



25 



Frequency response of ideal LMS TDL system. 

H3 m := E h t- exp (- j * fl, ^ T ) 
t 



30 



| H2 m | 0.5 



i 7 — r 




| H3 m | 0.5 



J L 




1 520 1 540 1 560 1 580 
X_Data m 



1520 1540 1560 1580 160 
X_Data m 
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10 



Gain compensation: 
Amplifier + TDL filter 
Comp k :=1(Mog[| A^-fl H2 |c |) 2 |] 



... Jc 

k 



15 



20 



25 



Normalize impulse response. 



Norm : 


t 


v= 


\ 




Norm 


|h« 1 


arg( h, ) 


0.654423 


0.014642 


0.05351 


2.334039 


0.135287 


-0.552664 


0.104525 


1.017118 


0.02811 


2.398415 


2.939684-1 O* 


2.781281 


0.01335 


2.126527 


7.851671-10"* 


-2992377 



24,169417 

max(Comp) - min(Comp) = 0.156938max Ripple in dB 



k 



E 1 

k 



= 0.047031 RMS Ripple 
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1 . An apparatus for changing the spectral profile of an optical signal comprising: 
a plurality of optical paths, each of said optical paths having associated 

5 therewith a delay element which imparts a predetermined delay to a 

signal propagating through the optical path; 
a plurality of couplers configured to split the optical signal among said plurality 

of optical paths according to predetermined splitting ratios; 
wherein said predetermined delays and predetermined splitting ratios are chosen 
10 so as to effect the desired change to said spectral profile. 

2. The apparatus of claim 1 wherein said predetermined delays and said 
predetermined splitting ratios are variable. 

15 3. The apparatus of claim 1 wherein said optical paths comprise waveguides and 

said delay elements include portions of said optical paths with an altered index of 
refraction. 

4. The apparatus of claim 3 wherein said altered index of refraction is permanently 
20 altered. 

5. The apparatus of claim 4 wherein said waveguides comprise silica waveguides, 
and said portions of said optical paths with an altered index of refraction comprise 
portions of said silica waveguides which have been exposed to UV radiation. 
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6. The apparatus of claim 2, wherein said waveguides comprise waveguides in 
electro-optic material and said delay elements include electrodes adjacent said 
waveguides. 

5 7. The apparatus of claim 3, wherein said waveguides comprise silica waveguides 

and said delay elements include heating elements adjacent said waveguides. 

8. The apparatus of claim 3 wherein said index of refraction is selectively 
reversibly alterable. 

10 

9. The apparatus of claim 3, wherein said optical paths comprise waveguides 
forming a splitter tree. 

10. The apparatus of claim 8 wherein said optical paths comprise waveguides 
1 5 forming a splitter tree. 

1 1 . The apparatus of claim 1 0 wherein said splitter tree has an input end and an 
output end, and further comprising a reflector at said output end. 

20 12. The apparatus of claim 9 wherein said splitter tree has an input end and an 

output end, and further comprising a reflector at said output end. 

13. The apparatus of claim 1 further comprising an optical amplifier and an optical 
path between said optical amplifier and said coupler. 

25 

14. A gain flattening filter for use with an optical amplifier comprising: 

an input waveguide; 
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a plurality of tapped delay lines comprising waveguides and delay 
elements associated with said waveguides, 

a plurality of variable couplers configured to split an optical signal 
propagating on said input waveguide among said plurality of 
tapped delay lines according to predetermined splitting ratios; 

wherein each of said delay elements is configurable to impart a 
predetermined delay to the signal propagating down, the 
associated waveguide, wherein said predetermined delays are 
chosen so as to effect a flattening of the spectral profile of said 
optical signal. 

1 5 . The gain-flattening filter of claim 1 4 wherein said tapped delay lines comprise 
waveguides formed in a planar substrate. 

16. The gain-flattening filter of claim 14 wherein said delay elements comprise 
selectively actuated resistive heating elements near said associated waveguides. 

17. The gain-flattening filter of claim 14 wherein said delay elements comprise 
regions of said associated waveguides having a propagation coefficient which differs 
from that of the remainder of said associated waveguide. 

18. The gain-flattening filter of claim 16 wherein said regions of said associated 
waveguides comprise regions which have been exposed to UV radiation. 

19. An optical communication system comprising: 

an optical amplifier; 

a spectral filter comprising an input waveguide; 
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a plurality of tapped delay lines, each of said tapped delay lines 
including a delay element on said line, 
a plurality of couplers configured to split an optical signal propagating 
on said input waveguide among said plurality of tapped delay 
5 lines according to predetermined splitting ratios; 

wherein said delay elements are configurable to impart a predetermined 
delay to the signal propagating down said line, wherein said 
predetermined delays are chosen so as to effect the desired 
change to said spectral profile; and 



10 an optical path between said amplifier and said filter. 

20. A method for modifying the spectral distribution of an optical signal comprising 
the steps of: 

splitting the signal among a plurality of optical paths, each having an 
15 output end, whereby a portion of the signal propagates on each of 



the optical paths; 

changing the relative phases of the portions of the signal propagating 

through the optical paths; 
modifying the amplitudes of the portions of the signal propagating 
20 * through the optical paths; • - - 

the relative phase changes and the amplitude modifications having 

magnitudes which result in a desired spectral distribution when 
the portions of the signal at the output ends of the optical paths 
are combined; and 

25 combining the portions of the signal at the output ends of the optical 

paths. 
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2 1 . The method of claim 20 wherein the step of changing the relative phases 
comprises delaying the propagation of the portion of the signal in at least one of the 
optical paths by an amount different from the others of said optical paths. 

22. The method of claim 21 wherein the amplitude modifying step comprises 
applying weighting coefficients to the portions of the signal propagating through the 
optical paths. 

23. The method of claim 22 wherein the delaying step further comprises modifying 
the refractive indices of each of said plurality of the optical paths. 

24. The method of claim 23 further comprising the steps of: 

determining the magnitude of each of said delays necessary to result in the 
desired spectral distribution; and modifying the refractive index of each of the 
optical paths by an amount which results in the determined delay magnitude. 

25. A method for modifying the spectral distribution of an optical signal comprising 
the steps of: 

splitting the signal among a plurality of optical paths having output ends 

according to predetermined splitting-ratios, 
changing the relative phases of the portions of the signal propagating through 

different ones of the optical paths; 
the relative phase changes and the splitting ratios having magnitudes which 

result in a desired spectral distribution when the portions of the signal at 

the output ends of the optical paths are combined; and 
combining the portions of the signal at the output ends of the optical paths. 



WO 02/14934 



PCT/US01/19788 



1/4 



FIG. 1 




38 



WO 02/14934 



PCT/US01/19788 



2/4 



FIG. 2 




WO 02/14934 



PCT/US01/19788 



3/4 



FIG. 4 






32-i 




31- 


CO 


30- 






o 


29- 


■s 




o 
5= 


28- 


I 


27- 




26- 




25- 



FIG. 5 




5 dB Ripple 



"1528 1533 1538 1543 1548 

Wavelength (nm) 



1553 



1558 



WO 02/14934 



PCT/US01/19788 




INTERNATIONAL SEARCH REPORT 



International application No. 
PCT/US01/19788 



A. CLASSIFICATION OF SUBJECT MATTER 
IPC(7) : G02F 1/01 

US CL : 385/1 

According to International Patent Classification (IPO or to both national classification and IPC 

B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 
U.S. : 385/1, 2, 3, 4, 5, 6, 7, 8, 9. 10, 14. 24, 47, 48, 49, 129, 130, 131, 132 



Documentation searched other than minimum documentation to the extent that such documents are included in the fields searched 
NONE 



Electronic data base consulted during the international search (name of data base and, where practicable, search terms used) 
Please See Continuation Sheet 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category ; 



Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



X 
X 
A 
X 
X 



US 4,840,447 A (KATAOKA) 20 June 1989 (20.06.1989), see entire document, especially 
Figure 1. 

US 5,253,309 A (NAZARATHY et al.) 12 October 1993 (12.10.1993), see entire 
document, especially Figure 9 and column 6, lines 26-28. 

US 5,347,601 A (ADE et al.) 13 September 1994 (13.09.1994), see entire document. 

US 5,351,317 A (WEBER) 27 September 1994 (27.09.1994), see entire document, 
especially column 2, line 62, through column 4; line 57, and Figures 1 and 2. 
US 5,832,155 A (RASCH et al.) 03 November 1998 (03.11.1998), see entire document, 
especially column 10, lines 3-38; column 14, line 59, through column 15, line 25; and 
Figure 8. 



1, 3, 8-10 
1-3. 6-10 

1-25 

1-25 

1-3, 6, 8-10, 13, 20-25 



□ 



Further documents are listed in the continuation of Box C. 



□ 



See patent family annex. 



* Special categories of cited documents: 

"A" document defining the general state of the art which is not considered to be 
of particular relevance 

"E" earlier application or patent published on or after the international tiling date 

document which may throw doubts on priority claim(s) or which is cited to 
establish the publication date of another citation or other special reason (as 
specified) 

u O" document referring to an oral disclosure, use, exhibition or other means 
"P" document published prior to the international filing date but later than the 



T" later document published after the international filing date or priority 

date and not in conflict with the application but cited to understand the 
principle or theory underlying the invention 

"X" document of particular relevance; the claimed invention cannot be 

considered novel or cannot be considered to involve an inventive step 
when the document is taken alone 

"Y" document of particular relevance; the claimed invention cannot be 

considered to involve an inventive step when the document is 
combined with one or more other such documents, such combination 
being obvious to a person skilled in the art 

**&" document member of the same patent family 



Date f the actual completion of the international search 
30 August 2001 (30.08.2001) 


Date of mailing of the international search report 

19 SEP 2001 


Name and mailing address of the ISA/US 
Commissioner of Patents and Trademarks 
Box PCT 

Washington, D.C. 20231 

Facsimile No. (703)305-3230 


"^•"S^ '.tiffin 

AkmUHah * 9 
Telephone No. 703-308-0956 



Form PCT/ISA/210 (second sheet) (July 1998) 



INTERNATIONAL SEARCH REPORT 



International application No. 
PCT/USOl/19788 



C ntinuation of B. FIELDS SEARCHED Item 3: 

USPTO EAST 

search terms: optical waveguides, splitting ratio, coupling ratio, couplers, phase, power, gain, flattening, equalize 



Form PCT/ISA/210 (second sheet) (July 1998) 



P.B.5818 - Patenllaan 2 
2280 HV Rijswijk <ZH) 
S +31 70 340 2040 
TX 31651 epo nl 
FAX +31 70 340 3016 



Europaisches 
Patentamt 

Zweigstelle 
in Den Haag 
Recherchen- 
ahtftiiunp 



European 
Pat nt Office 

Branch at 
The Hague 
Search 

division 



Office europeen 
des brevets 



Departemeni a 
La Haye 
Division de la 
iBchfirchfi 



r Brose, Gerhard 
Alcatel , 

Intellectual Property Depart ment Stuttgart/ 



70430 Stuttgart 
ALLEMAGNE 



Eingang bei ZPL 



17.FEa2003 



Datum/Date 

14.02.03 



L 


_l 1 1 


Zefchen/RefVRef. 

113 502 


Anmeldung Nr./Application NoiDemande n° ./Patent Nr ./Patent No./Brevel n". 

02360249.3-2415- 


melder/Applicant/Demandeur/Patentinhaber/ProprietoryTitulaire 

I ALCATEL 



COMMUNICATION 

The European Patent Office herewith transmits as an enclosure the European search report for the 
above-mentioned European patent application. 

If applicable, copies of the documents cited in the European search report are attached. 

□ Additional set(s) of copies of the documents cited in the European search report is (are) enclosed 
as well. 

The following specifications given by the applicant have been approved by the Search Division: 
El abstract QD title 

□ The abstract was modified by the Search Division and the definitive text is attached to this 
communication. 



The following figure will be published together with the abstract: 1 



REFUND OF THE SEARCH FEE 

If applicable under Article 10 Rules relating to fees, a separate communication 
from the Receiving Section on the refund of the search fee will be sent later. 




European Patent 
Office 



EUROPEAN SEARCH REPORT 



Application Number 

EP 02 36 0249 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Citation of document with indication, where appropriate, 
of relevant passages 



WO 02 14934 Al (WIGLEY PETER 6 ET AL) 
21 February 2002 (2002-02-21) 
* page 5, line 15 - page 6, line 15; claim 
1; figure 1 * 



US 5 132 639 A (FRAME JOHN S ET AL) 

21 July 1992 (1992-07-21) 

* column 2, line 42 - column 3, line 

* column 3, line 51 - column 4, line 

EP 1 204 228 A (CIT ALCATEL) 
8 May 2002 (2002-05-08) 

* abstract; figure 1 * 

US 5 812 294 A (WILSON GORDON COOK) 

22 September 1998 (1998-09-22) 

* abstract; figure 2 * 



1,7 



3 * 
39 i 



US 5 327 279 A (MCBRIEN GREGORY J 
5 July 1994 (1994-07-05) 
* abstract; figure 2 * 



ET AL) 



The present search report has been drawn up for all claims 



Relevant 
to claim 



CLASSIFICATION OF THE 
APPLICATION (lnLCI.7) 



H04B10/155 



1-12 
1-12 



1-12 



TECHNICAL FIELDS 
SEARCHED (lnt.CI.7) 



H04B 



Race o1 search 

MUNICH 



Date of completion of the search 

5 February 2003 



Phillips, S 



CATEGORY OF CITED DOCUMENTS 

X : particularly relevant if taken alone 

Y : particularly relevant if combined with another 

document of the same category 
A : technological background 
O : non-written disclosure 
P : intermediate document 



T : theory or principle underlying the invention 
E : earlier patent document, but published on, or 

after the Wing date 
D : document cited in the application 
L : document cited for other reasons 

& : member of the same patent family, corresponding 
document 



ANNEX TO THE EUROPEAN SEARCH REPORT 
ON EUROPEAN PATENT APPLICATION NO. 



EP 02 36 0249 



This annex lists the patent family members relating to the patent documents cited in the above-mentioned European search report. 
The members are as contained in the European Patent Office EDP file on 

The European Patent Office is in no way liable for these particulars which are merely given for the purpose of information. 

05-02-2003 



Patent document 
cited in search report 



Publication 
date 



Patent family 
member(s) 



Publication 
date 



wo 


0214934 26 


Al 








M AM C 

NUNL 








us 


5132639 


A 


21- 


07- 


Vive 


1 IC 


4992754 


A 


1 9—09—1 QQ1 

i.C \JC. 1331 














AU 


645878 


B2 


97— (11 —1 QQA. 














A 1 1 

AU 


1078292 


A 


n7— ni —i qq^ 














f* A 

LA 


2060117 


Al 


no— ns— i QQ9 














nr 

Ut 


69220514 


Dl 


^1—07—1007 
jl U/ 1:717/ 














ut 


69220514 


T2 


(1Q— 1 fl— 1 QQ7 
U7 1U 1:7:7/ 














tr 


0498456 


Al 


1 9— Oft— 1 QQ9 

it UO l73t 














Lb 


2104741 


T3 


1 (\— 1 fl— 1 Q07 
ID 1U 137/ 














JP 


3221594 


B2 


22-10-2001 














JP 


7202763 


A 


04-08-1995 














ID 

or 


2056370 


C 
















JP 


4336819 


A 


25-11-1992 














JP 


7087403 


B 


20-09-1995 














us 


5363056 


A 


08-11-1994 














US 


5252930 


A 


12-10-1993 














A 1 1 

AU 


645029 


B2 


Hfc— 01 — 1 QQA 
UO Ul 117 74 














A 1 1 

AU 


6199590 A 


1 A fi^— 1 QQ1 

m— uj- iyy i 














n a 
CA 


2024385 


Al 


nQ-n'i-i QQi 
uo— uo iyy i 














DE 


69032831 


Dl 


9Q_m _1 QQQ 

do— ui— 














Ut 


69032831 


T2 


1 o_ QQQ 

XL. UO 1^7i7i7 














DE 


416622 


Tl 


19-03-1992 














EP 


0416622 


A2 


13-03-1991 














ES 


2125851 


T3 


16-03-1999 














JP 


1925165 


C 


25-04-1995 














JP 


3179807 


A 


05-08-1991 














JP 


6052816 


B 


06-07-1994 


EP 


1204228 


A 


08- 


-05- 


-2002 


EP 


1204228 


Al 


08-05-2002 














US 


2002054421 


Al 


09-05-2002 


US 


5812294 


A 


22 


-09- 


-1998 


NONE 








us 


5327279 


A 


05 


-07 


-1994 


CA 


2144960 


Al 


03-02-1994 














DE 


69316867 


Dl 


12-03-1998 














DE 


69316867 


T2 


28-05-1998 














EP 


0612452 


Al 


31-08-1994 














WO 


9402998 


Al 


03-02-1994 




o 

Si For more details about this annex : see Official Journal of the European Patent Office, No. 12/82 



